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We analyze and clarify the transport properties of a one-dimensional metallic nanoparticle array 
with interaction between charges restricted to charges placed in the same conductor. We study the 
threshold voltage, the I-V curves and the potential drop through the array and their dependence on 
the array parameters including the effect of charge and resistance disorder. We show that very close 
to threshold the current depends linearly on voltage with a slope independent on the array size. 
At intermediate bias voltages, for which a Coulomb staircase is observed we find that the average 
potential drop through the array oscillates with position. At higher voltages I-V curves are linear 
but have a finite offset voltage. We show that the slope is given by the inverse of the resistances 
added in series and estimate the voltage at which this linear regime is reached. We also calculate 
the offset voltage and relate it to the potential drop through the array. 
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Nanoparticle arrays made of metallioii^i^ii^i^i^'^'^'i^, 

seniiconductin g^°i^^d^i^^d"^i^^i^^ , magneti o^^i^^'^^ or 
combined2£i2iiS^ materials and with radii of the order 
of 2-7 nm can be now synthesized. The transport 
properties of these systems are influenced by the 
ratio between the energy level spacing, the charging 
energy of the nanoparticles, and the temperature. 
The first two quantities depend on the material and 
the size of the nanoparticle. In the case of metallic 
nanoparticles, at not too low temperatures, the level 
spacing is much smaller than the temperature and 
does not play any role in the transport'^'^. On the 
contrary, the charging energy is of the order of 0.1 eV. 
Strong interactions between the electric charges and 
the possibility of tuning interparticle coupling make 
nanoparticles arrays an ideal system to study correlated 
mot ion^'^ -^^-^^-^'^-^^-^Q-'^O-^^-^^''^^-^^'^^-^^''^''-^^-'^^-^"-^^-'^^ 

Experimentally, these arrays are strongly influenced 
by disorder— i^i^. Local charging disorder is present 
in all arrays due to randomly dispersed charged impu- 
rities lodged in the substrate or in the materials that 
separate and surround the nanoparticles. Because of the 
exponential dependence of the tunneling resistance, even 
a small dispersion in the distance between nanoparticles 
results in large variations in the tunneling resistances of 
the junctions. Differences in the island sizes and voids in 
the lattice can be other sources of disorder—. 

Due to the combination of disorder and charging ef- 
fects the current in voltage biased arrays is blocked up to 
a threshold voltage^i-^'i^ii&iiiii^iiai^'^'^ Vt- For bias 
voltages larger than Vt current is in general non-linear 
in voltage with a power-law dependence^i^iS close to 
threshold, a linear dependence recovered at high- voltages 
and frequently a step-like behavior, called a Coulomb 
staircase, at intermediate voltages. Most studies have 
focused on the statistical analysis of the threshold volt- 
age and on the power-law behavior of the current close 



to this threshold. This exponent depends on the dimen- 
sionality of the array, but there is controversy between 
different theoretical approaches in the one-dimensional 
case, with both linea r^^i^^i^^ and square-root^ predic- 
tions. Much experimental work has been concentrated 
on two and three dimensional arrays, but some quasi- 
one dimensional systems have also been fabricate d^'^'^^ . 
Comparison between experiments and theory is not yet 
well settled. 

In this paper we provide a complete description 
of the zero-temperature transport properties of one- 
dimensional metallic nanoparticle arrays for the case in 
which interactions are restricted to charges in the same 
nanoparticle (onsite limit). In a separate work^ we dis- 
cuss the effect of the long-range character of the inter- 
actions on these properties. We resolve the controversy 
on the power-law, analyze carefully dependencies on dif- 
ferent array parameters and estimate the voltage which 
delimits each of the transport regimes to favor compari- 
son with experiment. We discuss arrays with and with- 
out charge disorder. Although clean arrays are mainly 
of academic interest, their analysis will help us to un- 
derstand the main features of the experimentally more 
relevant, disordered arrays. The effect of variations in 
the junction resistances is also analyzed. Disorder in 
capacitances (nanoparticle size variations) is not con- 
sidered as it is less important in present experiments. 
Nanoparticles synthesized nowadays are monodispersed 
in size to a few percent. In any case, the effect of capac- 
itance disorder— in most of the properties studied here 
can be deduced from the analytic approximations pro- 
vided in the text. Due to the one-dimensionality of the 
array and the nearest neighbor tunneling considered we 
assume that there are no nanoparticle voids in the array, 
as this would completely prevent current flow. We have 
analyzed the threshold voltage, the TV characteristics, 
close to threshold but also at larger bias voltages and 
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the potential drop through the array. Among the results 
presented we show that the threshold voltage of clean ar- 
rays in the onsite case does not approach the one found 
for capacitively coupled junctions in the limit of weak 
coupling^i^. The controversy in the power-law depen- 
dence of current on bias voltage close to threshold is clar- 
ified and shown to be linear, but with a slope different 
from the one previously predicted§4. We find signatures 
of correlated charging of the system in the potential drop 
through the array, especially at intermediate values of the 
bias voltage. We show that even at high voltages, with 
the conductance controlled by the sum of the resistances 
in series, the potential drop through the array is not com- 
pletely linear, but it has a contribution due to charging 
effects. We identify an asymmetry external parameter a 
which controls the bias voltage drop. The influence of a 
has barely been discussed in previous works. We show 
that its effect can a priori be observed experimentally in 
the width of the Coulomb staircase steps. Calculations 
are perfomed numerically, but analytic approximations 
are given in several limits and compared with numerical 
results. 



I. THE MODEL 

We analyze the current through a one-dimensional 
array of N metallic nanoparticles placed between two 
electrodes. We consider the classical Coulomb blockade 
regime with S ^ KbT < El'^K S is the level spacing, T 
the temperature and = 1/(2C*^') the charging en- 
ergy of the islands, with capacitance C^^'. Here and in 
the following the electronic charge e = 1 We assume that 
each nanoparticle has a continuum level spectrum (S — 0) 
and a constant density of states at the Fermi level, but 
a gap E"^^^ for adding charge. We restrict electrostatic 
interactions to those charges on the same conductor: ca- 
pacitive coupling vanishes. This limit is referred to as 
short-range or onsite interaction limit. The nanoparti- 
cles are separated by high tunneling barriers with a re- 
sistance much larger than the quantum of resistance. In 
these conditions the charge in the islands can be assumed 
fixed and quantized. Eventually we allow tunneling pro- 
cesses between nearest neighbors, and treat the transport 
at the sequential tunneling level. A single charge is in- 
volved in the tunneling process. We assume that when a 
charge hops, the charge density in the final state of the 
array immediately relaxes to the electrostatic equilibrium 
configuration. 

We take into account that the electrodes are not ideal 
voltage sources, but have a finite self-capacitance. In 
equilibrium, and before the tunneling event the elec- 
trodes are held at a given potential due to the charge 
provided by a battery. We assume that the tunneling 
time, i.e. the time needed by the electron to cross the 
tunnel barrier, is smaller than the circuit characteristic 
time that determines how quickly the battery can trans- 
fer charge to the leads in order to restore the voltage at 



the electrodes. As a consequence, just after the tunneling 
process the electrodes will not necessarily be at the same 
potential at which they were at before the tunneling event 
because the charge, provided by the battery, necessary to 
restore their initial potentials has not arrived yet. The 
voltage is restored to the nominal value before the next 
tunneling event. For finite-range interactions the poten- 
tials on the leads will thus fluctuate in response to all 
tunneling events, even those that do not directly involve 
the electrodes. In the short-range case, considered here, 
they will fluctuate only when an electron jumps into or 
out of the leads. 

The current is calculated numerically by means of 
a Monte Carlo simulation, described in the Appendix, 
which depends on the tunneling rates. The probability 
of a tunneling process'^^ is given by 



r(A£;) 



1 



R exv{l\ElKBT) - 1 



(1) 



with i?, the tunneling resistance of the junction. We 
will restrict the discussion to zero temperature for which 
r(A£;) = -A£:/i?e(-A£;). as is the difference be- 
tween the energy of the system before and after the tun- 
neling event, with the sign convention that Ai? is nega- 
tive if the energy decreases. It excludes the work done 
by the battery to recharge the electrodes, as explained 
before. The energy gained by tunneling is assumed to 
be dissipated. Only changes in energy with electrostatic 
origin are considered. The energy of our system is given 
by 
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Labels and TV -|- 1 refer to source and drain electrodes 
and 1,...,A^ to the islands. In the following, latin capital 
and lower case letters are used to denote electrodes and 
islands respectively. Greek indexes will be used when the 
labels refer to both islands and electrodes. Ca and Qa are 
respectively the conductor capacitance and the charge in 
it. The charges provided by the battery at the source 
and drain electrodes, which maintain their potentials at 
Vo and Va^+i are Qq = CqVq and Qn+i = Cn+iVn+i- 
The capacitances of the electrodes are much larger than 
those of the nanoparticles. c/if'" is a random potential 
at each island due to randomly dispersed charges within 
the substrate and within the material surrounding the 
nanoparticles^. Clean arrays will be characterized by 
^dis _ Q gvery i. The random potential is included 
only at the islands because a similar term at each elec- 
trode is compensated by the battery and thus has no 
effect on transport. In the case of disordered arrays, the 
disorder potentials can, in principle, take values larger 
than the charging energy E^^K However, for large val- 
ues of the disorder potential, charges flow to compensate 
for these large fluctuations. In the case of short-range 
interactions, except if the original disorder potential is 
very weak, once the screening of the potential due to 
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the mobile charges is taken into account, the set of dis- 
order potentials is uniformly distributed in the interval 
__^j,s; ^dis ^ ]t;^s143^ g^jjjj ^Yie following we consider 

this distribution. 

The relevant quantity for the transport is the change 
in energy due to a tunneling event. The tunneling pro- 
cess can be seen as the creation of a hole in the conductor 
a from which the charge leaves, Qa — > Qa — 1, and the 
addition of an electron in (3 at which the charge arrives, 
Q/3 — * Q/? + 1- Here and thereafter, we let 1) de- 

note the charge of an electron (hole). In fact, the change 
in energy can be rewritten as the energy to create an 
electron-hole (also called in the following excitonic en- 
ergy) plus the difference in potential between the sites 
involved in the process before the tunneling event. 



AE = E'^J + {ct^p - cj,^) 



(3) 



The first term gives the energy to create an electron-hole 
pair in an uncharged clean array and is given by 
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This energy does not depend on the direction of tunneling 
(from a to /9 or from /3 to a) and in the following it will 
be denoted E^~^ with i running from 1 to + 1 . Index i , 
when used to label a junction will refer to the one between 
conductors i — 1 and i. We will use the term contact 
junction for those junctions which connect an island and 
an electrode, and bulk or inner junction for those ones 
in between two nanoparticles. For the contact junctions 



i = \,N + ie: 



e — h 



E'f' as E. 



source, drain 



« El"^ while 



for the bulk junctions i^2toi^ N E^^'' = 2El''K The 
second term in ([3]) can be seen as the change in potential 
between the conductors involved in the tunneling. The 
potential at each site depends on the charge state of the 
array prior to the tunneling event. At the electrodes (j)o = 
Vb, (t>N+i — Vn+i- At the islands, the potential can be 
decomposed into two terms: (f>f'^'' , a random potential due 
to random charges in the substrate and (p^^ a potential 
due to the charges in the islands. ~ (j)f^^ + (jjf^ , with 
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Analogously we can define the potential drop at each 
junction 



(6) 



with the corresponding disorder and charging terms ^f"^^ 
and <I>^'\ The potential drop at a contact junction de- 
pends on the disorder and charge state of the nanoparti- 
cle and on the applied bias voltage. On the contrary the 
potential drop at a bulk junction is not affected by the 
bias voltage, except via a change in the charge state. 

In the following we rewrite the potential at the elec- 
trodes as Vb = aV and V^+i = {a — 1)V. The total 
potential drop through the array is Vq — V^+i = V . In 



our model some measurable properties depend on the 
value of a, which characterizes how the bias voltage is 
partitioned between source and drain chemical potential 
shifts. In several previous works the value of a was cho- 
sen cither as a — 1/2, correspondingly to a symmetrically 
biased array, or as a = 0, 1 corresponding to completely 
asymmetric biasing, a — 1 has been also called the for- 
ward bias condition^^ Both values have been used in the 
literature, mostly without discussion. In the symmet- 
rically biased case the potential drop at both contact 
junctions is equally modified by the bias voltage. On the 
contrary for a — 0(1) only the drain (source) junction is 
affected by the bias. Since no physical properties depend 
on the overall zero of energy, varying a in our model is 
entirely equivalent to rigidly shifting all impurity poten- 
tials by —aV. The dependence on a discussed below, 
corresponds in part to a dependence on the alignment of 
the equilibrium source and drain chemical potentials with 
respect to the addition and removal energies of the elec- 
tron. For example for a single nanoparticle on whether 
the chemical potential shift required to add or remove an 
electron is larger. Since in our model all transport occurs 
by transfer between adjacent nanoparticles, the evolution 
of a nanoparticle array as the bias voltage is applied is 
sensitive to a. For a given nanoparticle array with a fixed 
set of disorder potentials, we believe that the dependence 
on a discussed below should in principle be observable. 

Whenever not specified we assume that all the junction 
resistances Ri are equal and given by Rt- The effect 
of non homogeneous resistances will be studied in two 
ways. One of the junction resistances at a given position 
is larger than the other ones (given by Rt) or resistances, 
varying in between two values are randomly assigned to 
the junctions. To mimic that disorder in resistances orig- 
inates in variations in distances between the islands and 
the exponential dependence of the junction resistance on 
the distance between islands the junction resistance is 
given by i? = RQexp{"fdist) with Rq and 7 input pa- 
rameters and dist = 1 -t- random/ 2. Here random is a 
random number between and 1. In the paper, we have 
used Ro = 1.1825i?T and 7 = 1.526,1.95,2.84. With 
these values the resistance changes respectively between 
(5-ll)i?T, (8-21)i?T and (23-83)i?T. 



II. THRESHOLD VOLTAGE 

In this section we analyze the dependence of the 
threshold voltage on the array parameters. The thresh- 
old voltage is controlled by changes in energy in tunneling 
and not affected by the resistance of the junctions. Thus, 
we do not address the case of disorder in resistance in this 
section, as it does not modify the threshold compared to 
the equal-resistances case. 

The threshold voltage is the minimum bias voltage at 
which current can flow through the array. A finite bias 
voltage can assist the entrance of charge to the array 
from the leads, as it creates a potential drop at the con- 
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tact junction, which can overcome the excitonic energy. 
A finite current requires that charges are able to be trans- 
ferred from one electrode to the other one across the en- 
tire array. If charge flow can occur between the leads, the 
threshold voltage is the minimum voltage which permits 
the entrance of an electron or hole into the array. How- 
ever it is possible for charge to become stacked inside the 
array due to the disorder potential configuration or due 
to the lack of potential drops across the bulk junctions. 
In this case the threshold will be controlled by the flow 
of charges. 

Here for disordered arrays we recover previously 
predicted^^ values for the average threshold < Vr > and 
its root mean square deviation 6Vt- < Vt > is propor- 
tional to the number of particles and SVr to N^^'^ . On the 
contrary, for clean arrays and strictly onsite interactions 
the threshold voltage differs from the one expected by ex- 
trapolating to zero coupling the value obtained for weakly 
coupled nanoparticles. In particular, for finite coupling 
(or its zero coupling extrapolation) an A^-independent 
threshold voltage is predicted for large arrays, while for 
strictly zero coupling we obtain Vr ~ 2Eq''N. We also 
show that its value depends on a. 

The threshold corresponding to the clean case is plot- 
ted in Fig. 1(a) for the case of symmetrically biased ar- 
rays (a = 1/2), antisymmetrically biased arrays {a — 1), 
and an intermediate biasing (a = 3/4). In the symmet- 
rically biased case, Vt shows a step-like dependence on 
N. There is a clear even-odd effect. The even-odd effect 
in the threshold voltage with the number of particles is 
absent for a = 1 and a = with a threshold voltage 
Vt = El^'-i2N — 1). At strictly zero temperature the 
tunnehng rate T{/S.E) — ■^AEQ{—AE) and vanishes 
when AE is positive or zero. At the inner junctions AE 
is independent of the bias voltage and zero or negative if 
the two islands differ just by just a single charge. Thus 
the tunneling rate vanishes. A^ — 1 junctions prevent the 
flow of charges. A charge gradient at each bulk junction 
has to be created to allow flow of charge. In the symmet- 
rically biased case, increasing the potential at the elec- 
trodes allows positive and negative charges to enter from 
the source and the drain respectively. These charges ac- 
cumulate on the array and create potential drops across 
the bulk junctions. At voltages just below the threshold, 
the accumulated charges at the first and last islands are 
equal in number and opposite in sign. Current starts to 
flow at voltages larger than Vr — 2NEl^^ when A^ is odd 
and Vt = 2{N — 1)E].^^ when A^ is even, as these values 
allow for the build up of ±{N - l)/2 and ±{N - 2)/2 
charges at the first and last islands for odd and even 
N respectively. This corresponds to a charge gradient 
dQi — Qi — Qi-i = —1, across all bulk junctions for odd 
A^ and across all bulk junctions except one for even A^. 
For a = 0, 1, charge can only enter the array from one 
lead and the energy barriers across all N — 1 bulk junc- 
tions must be overcome by accumulated charges. When 
a = 1/2 the first and last junctions are equivalent and 
charge can enter from both leads so it is possible in some 



cases for the energy barrier across one of the bulk junc- 
tions to be overcome by the potential drop due to two 
injected charges of opposite sign from the two leads, i.e. 
one of the junctions can be uncharged. The absence of 
this possibility is what removes the even-odd effect when 
a = and a = 1. An intermediate situation is found for 
a = 3/4. 

Complementary information can be obtained by look- 
ing at the threshold voltage as a function of a for different 
values of A^ in Fig. 1(b). The threshold voltage changes 
in a periodic way with a. The periods of the features 
in Vt (a) depend on the number of barriers in the array. 
Dependence of Vt on N and periodic features in Vr with 
respect to a reflect the number of charges which have to 
accumulate in the flrst and last island prior to current 
flow. Fig. 1(b) includes curves corresponding to Vq and 
Vn+1 for speciflc values of Qi and Qn. As the assym- 
metry of the array increases (increasing \a — 1/2|), the 
threshold voltage is alternately determined by the cost 
of injecting a charge unto the array from the source and 
the drain. At all a except for the values that lead to the 
minimum Vt for a given array length, the difference in 
the charge occupying the first and last islands, Qi~-Qn, 
equals A" — 1. At the minimum values of Vr, this charge 
difference equals N — 2. 

The clean case of a related system was studied by Hu 
and O'ConneflSS. They analyzed a one-dimensional array 
of N gated junctions with equal junction capacitances Cj 
and equal gate capacitances Cg. Due to the finite value 
of Cj charges in a given island interact with charges in 
other islands and with charges in the electrodes. With 
an applied bias voltage the interaction between charges in 
the electrodes and in the islands results in a bias induced 
potential drop at the bulk junctions. Once a charge is 
injected unto the array, it will have no difficulty in trav- 
eling through it, and the threshold voltage equals the 
voltage required for injection of a charge from the elec- 
trodes. As the ratio Cg/Cj increases the threshold volt- 
age of a long array tends to an N-independent value of 
the order of the charging energy. The onsite case dis- 
cussed here corresponds to Cj = 0. If one extrapolates 
the case discussed by Hu and O'Connell^ to Cg/Cj — > 
an N-independent threshold voltage would be expected 
for onsite interactions. As shown above, the threshold 
voltage of clean arrays does not satifies this prediction, 
as at zero temperature the charges cannot travel freely 
through the array, and the threshold voltage increases 
with the number of islands. 

In the case of disordered arrays, each array has a 
threshold voltage Vr dependent on the given configura- 
tion of disorder {0''**}. The threshold voltage depends on 
a, in a way similar to the clean case, see Fig. 1(b). How- 
ever, as shown in Fig. 1(c), this dependence on a disap- 
pears in the average value and we recover Middleton and 
Wingreen prediction^. For the disordered case Middle- 
ton and Wingreen^l predicted a linear dependence of the 
threshold voltage on the array length. Upward steps in 
the disorder potential $f'* > prevent the flow of charge. 



5 



20 
16 

12 



~i — I — I — I — I — I — I — I — I — i7~r 

□ — H a=l , 




Q l . I I I I I I I I I il I I I I I I I I 

2 4 6 8 10 0.2 0.4 0.6 0.8 1 

N cx 

60 



40 

> 

V 20 



to 



1 



□ a=l 
* a=0.5 



T 



^ (c) 



10 




-- v^(n)=e;''n 

i'' □ a=l 

A a=0.5 
a=0.75 

\ ^ \ ^ 



20 40 
N 



60 



FIG. 1: (a) Threshold voltage of clean arrays as a function of 
the number of islands TV for different values of the asymme- 
try bias parameter a. A clear even-odd effect is present for 
symmetrically biased arrays, a = 1/2, (Vb = V/2 = — Viv+i). 
(b)Threshold voltage as a function of a. From top to bottom 
lines with symbols correspond to clean arrays with 6 and 5 is- 
lands and to disordered arrays with 6 and 5 islands. The peak 
and valley dependence reflect the change of contact junction 
which control the entrance and the number of charges which 
have to accumulate in the corresponding contact island. For 
comparison thin solid and dashed lines give, for clean arrays, 
the bias voltage at which 2 and -2 charges can be placed at 
the first and last islands, (c) Main figure: average threshold 
voltage for disordered arrays as a function of the number of 
islands. The dependence of Vt on a disappears on average 
and a linear dependence on A'^ is recovered. Fluctuations in 
threshold voltage follow SVt ~ N^^^ as predicted^ 



The downward steps < facilitate it. In average 
there are N/2 upward steps. To overcome such steps a 
charge gradient has to be created in those junctions. For 
onsite interactions this results in^ <Vt>— E^^^N . In 
this limit, they also argued that fluctuations in the values 
of Vt over many random configurations of disorder are 
analogous to the fluctuations in the net distance trav- 
eled by a 1-D random walk with N steps, resulting in 

Compared to arrays with no disorder, the threshold of 
disordered arrays in the onsite limit tends to be smaller 



because in the clean case, all the bulk junctions inhibit 
charge flow whereas in the disordered case, only the bulk 
junctions with positive <l>f** block charge flow. As shown 
in the inset of Fig. 1(c), we also recover the relation- 
ship for the fluctuations in the threshold predicted by 
Middleton and Wingreei*^, 5Vt oc N^/'^. 



III. FLOW OF CURRENT 

For bias voltages larger than threshold the current / 
can flow, but it is a strongly non-linear function of volt- 
age. The current depends on the charging energy and 
number of islands, the presence or not of charge disorder 
in the array, the resistances of the junctions and on the 
asymmetry of the applied bias voltage. Linearity and 
independence on a is recovered at large voltages. As 
discussed by Middleton and Wingreeen^ a power-law 
dependence of the current on voltage is found close to 
threshold. In this section we discuss the different regimes 
which can be differentiated in an 1-V characteristics and 
its dependence on the input parameters. At very low 
voltages by comparing analytical and numerical results 
we resolve the controversy on the exponent of the power- 
law and show it to be linear with a slope which depends 
on a and the resistance of the contact junctions, but that 
is independent on the array length (except in a particular 
case in which an even-odd effect is found). The linearity 
is however restricted to very small values oiV — Vt- We 
also clarify the dependence of the Coulomb staircase pro- 
file on the bias parameter a and estimate the asymptotic 
current at high voltages, and the bias voltage at which 
this high voltage regime is found. Very large values of 
the bias voltage have to be applied to reach this linear 
dependence. 



A. Linear dependence close to threshold 

There has been some controversy regarding the power- 
law of the current with {V — Vt) through one dimensional 
disordered arrays for voltages close to Vt, / ~ {V — VtY ■ 
Middleton and Wingreenj^i predicted linear behavior for 
both the long and short range interaction. Reichardt and 
Reichardt^^ found a square root behavior using a model 
with a l/r interaction between the charges in the islands. 
They argued that = 1/2 is the exponent corresponding 
to an sliding charge-density wave. They pointed out that 
the larger values of the exponent obtained by Middle- 
ton and Wingreer^ are a consequence of using voltages 
which are not small enough. Kaplan et al.'^^ found // = 1 
in the long-range limit of an array of dots capacitivcly 
coupled to their nearest neighbors. Finally Jha^^ and 
Middleton argued that the dependence of the current of 
disordered arrays in the onsite limit on {V — Vt) for volt- 
ages marginally greater than Vt is linear with an slope 
inversely proportional to the length of the array. Numer- 
ically they found an approximate linear behavior only in 
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the case of very long arrays but not for the smallest volt- 
ages analyzed (where they found a sublinear dependence) 
but in an intermediate voltage regime. So far, there are 
no experiments available in completely one-dimensional 
arrays, but there are a few in quasi-one dimensional sys- 
tems. The approximate power-law measured'*- at volt- 
ages {V — Vt) ^ O.OIVt is larger than unity which has 
been attributed to the fact that the system is not strictly 
one-dimensional. 

In this subsection, we show that the current varies lin- 
early with respect to V — Vr for very small V — Vt and 
that the slope is not inversely proportional to N. The 
slope is sensitive to the resistance of the contact junc- 
tions and to the degree of symmetry between the applied 
bias voltages on the source and drain leads. While this 
linear dependence is found analytically and in the nu- 
merical calculations the value of {V — Vr) at which it 
disappears is very small. Because of the small value of 
the current at such voltages, the linear dependence most 
probably cannot be seen experimentally. 

The linearity above but very close to threshold can be 
easily understood. Charges can enter only through the 
contact junctions from the leads. The current through 
the array is equal to the average charge transferred per 
unit time. The average time necessary to transfer a 
charge through the array is the sum of the time in- 
volved in all the processes in the sequence of tunneling 
events from the moment in which charge enters the ar- 
ray from one electrode until it leaves the array to the 
other one. If the time associated to tunneling at a given 
junction is much larger than the time involved in the 
rest of processes, this junction acts as bottleneck and the 
time necessary to transverse the array is approximately 
equal to the inverse of the scattering rate through this 
junction. Thus, the current can be approximated by the 
tunneling rate across the bottle-neck junction. Below 
threshold, but close to it, the tunneling at the contact 
junctions costs finite energy and transport is suppressed. 
This cost in energy is reduced by the applied bias volt- 
age and at threshold is zero at the entrance junction. 
When V > Vt but very close to it, the junction from 
which the charge enters the array acts as bottle-neck for 
transport, once it is allowed. At zero temperature, the 
tunneling rate is —AE/R. The dependence on bias volt- 
age of the energy for tunneling compared to the one at 
threshold, at the source and drain junction is, respec- 
tively JiJs = -a{V-VT) and SEd = -{1 - a){V ~ Vt). 
If charge enters through a single junction, the current 
approximately equals 



I = -^aiV ~ Vt) 
Hi 

if charge enters from the source, and 
1 



Rn 



+1 



-{l-a){V -Vt) 



(7) 



(8) 



if holes enter unto the array from the drain. Both source 
and drain junctions have to be taken into account in the 
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FIG. 2: Main figures: I-V curves in logarithmic scales for 
different array parameters. Insets show the derivatives (in 
units of 1/Rt) of the curves plotted in the main figures. The 
linear dependence of current on voltage (constancy of the 
derivative) is clearly seen in all the plots but it disappears 
for V-Vt 10-'^ Ei"' which for the cases shown corresponds 
to {V — Vt)/Vt ~ 10~*. (a) to (c) show I-V curves cor- 
responding to arrays without charge disorder. All junction 
resistances are equal in (a) and (b). As shown in (a) for clean 
arrays the slope of the linear dependence does not depend 
on the number of islands except for a = 0.5 which shows an 
even-odd effect with a slope equal to unity and 0.5 for odd and 
even number of particles, respectively, (b) The dependence of 
slope on a can be non-monotonous if there is a change in the 
contact junction which acts as a bottle-neck, (c) I-V curves 
for arrays with junction resistances randomly assigned vary- 
ing between (5 — 11)7?t (upper curve), (8 — 21)Rt (middle 
curve) and (23 — 83) Rr (bottom curve). The different slopes 
are due to different resistances at the bottle-neck junction, 
(d) I-V curves of disordered arrays with homogeneous con- 
tact resistance. The a = 0.5, N = 5 curves correspond to 
different realizations of disorder. The even-odd effect present 
for a — 0.5 in the clean case has disappeared as just one 
contact junction acts as bottle-neck. 



clean symmetrically biased case when N is odd. 



Rn+i) 2^^ 



Vt) 



(9) 



When A'' is even and a = 0.5 it is necessary that charge 
enters through both junctions for current to flow and 
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current is approximately equal to 

i = irTj^ — \{y-yT) m 

This behavior is observed in Fig. 2. The linear behavior is 
clearly appreciated in both the log-log scale in which the 
main figures are plotted as well as in the constancy of the 
derivatives in the insets. The dependence of the slope of 
the I-V curves is better seen in the conductance plotted 
in the insets. As seen in (a) the dl/dV is equal for the 
a = 0.5, = 50 and = 20 curves, i.e. independent on 
the array length. On the contrary, it is double for A^ = 
5. This behavior originates in the even-odd alternancy 
predicted by Eqs.® and (jTU]). The dependence on a is 
studied in Fig. 2(b). From top to bottom, the change of 
slope with a is smooth and given by a/ Ri as predicted 
by ([7|), but turns non monotonously if currents starts 
being controlled by ([5]) with slope (1 — a) / Rn+i- As also 
expected from the discussion above the slope is affected 
by disorder in resistances in Fig. 2(c), but not by charge 
disorder in Fig. 2(d), except in the lose of the even-odd 
effect present in the clean case for a = 0.5. 

B. Loss of linear behavior and intermediate regime 

The linear behavior in Fig. 2 appears for several or- 
ders in magnitude. In spite of this, it disappears for 
V -Vt^ 10-2 or {V - Vt)/Vt ~ lO-"*. The magnitude 
of the current is probably too small for this linearity to 
be detected experimentally. We have found that linear- 
ity gives rise to sublinear behavior when the time of the 
other tunneling processes become relevant compared to 
the time spent at the bottle-neck. Upon increasing V the 
tunneling rates of the different processes involved in the 
transport become more homogeneous. To obtain sublin- 
ear behavior, it is just necessary than the two slowest 
processes in a sequence have comparable rates. The lost 
of linearity can depend on the resistance of the junctions 
when a non bottle-neck junction has a resistance much 
larger than the bottle-neck one. The bottle-neck charac- 
ter of a process at a junction at the leads can disappear 
faster for longer arrays as there are more tunneling pro- 
cesses which will contribute to the total time. In the 
disordered case, the energy gain of some of the tunneling 
processes is smaller than in the clean case and the con- 
tact junctions can stop being the bottle-neck earlier, i.e. 
for smaller V — Vr- But, in general we have not found 
very significative differences for different array parame- 
ters in the value of the bias voltage at which the linearity 
disappears. 

As the energy of a tunneling process through an inner 
junction does not depend on the applied voltage (except 
via the charge accumulated on it), when a bulk junc- 
tion controls the transport, the current is independent 
on voltage showing a characteristic staircase profile. The 
existence of the Coulomb staircase has been known for 
a long time^^ii^i^i^. Early claims reported a Coulomb 



staircase only in the asymmetrically biased case^. More 
recent results in clean capacitively coupled nanoparticle 
arrays, show that a staircase also emerges in a symmet- 
ric array under symmetric bias^^, but claim that the I-V 
characteristic for an N-dot array under forward bias is 
identical to that for a 2N-dot one under symmetric bias. 
We show here that while the appearance of the staircase 
is generic, the last statement is not correct. 

The current has kinks at those voltage values which 
change the maximum number of charges which can be ac- 
cumulated at the first or last island, allowing new trans- 
port processes. These voltage values depend on the asym- 
metry a and the existence of charge disorder in the array 
but not on resistance disorder. To allow the addition of 
an extra charge in the first or last nanoparticle requires 
an increase in bias voltage in the adjacent electrode of 
approximately 2i?**'. In the clean case, when a = 1,0 
only one electrode changes its potential and the width 
of the steps in bias voltage is 2i?*''' . On the other hand, 
when a = 1/2 the change in potential of a given elec- 
trode is just the half of the bias voltage and steps appear 
in intervals of 4i?*'*' . 

With charge disorder and a = 0, 1 the position of the 
kinks slightly depends on voltage, but the width of the 
voltage intervals between the kinks does not change, as 
new charges are added through a single junction. If a = 
1/2, charges enter from both contact junctions but the 
corresponding kinks in the current do not appear at the 
same position. While the width of a kink corresponding 
to a given junction remains 4i?**' in a general case in the 
I-V characteristic there will be two kinks in each 4i?*'*' 
interval in bias voltage due to the alternative position 
of the kinks of both contact junctions. Except in very 
special cases the separation of kinks does not do equal 

The position of the kinks in the clean and disorder 
cases can be observed in Fig. 3. Two main features can 
be observed. For clean arrays (insets in (a) and (b)) 
at the onset the current shows a big jump. Once the 
charge gradient is created and charge can transverse the 
array, it fiows easily. The steps at higher voltages are 
much smaller but have the width predicted above. In 
the main figures corresponding to disordered arrays, the 
large big jump has disappeared and steps are more clearly 
observed. Two stepwidths, which add AE'^'' are seen for 
a = 0.5 and just steps with width 2E^^ appear for a = 1. 

Previously, in the analysis of the threshold voltage, we 
saw that while the value of a has an effect the value for 
the onset of current, it does not seem possible to deter- 
mine a from a measurement of Vr in the, experimentally 
relevant case, of a disordered array. However if the on- 
site interaction case discussed here can be experimentally 
reproduced, the width of the steps in the I-V curve can 
differentiate a values. Disorder in the resistances does 
not modify the voltages at which kinks in the current 
appear but it does affect the staircase profile. The stair- 
case profile is modified in Fig. 2(c), compared to (a) due 
to the disorder in resistances. A very large resistance 
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FIG. 3: I-V curves for N=50 and different array parameters 
at intermediate bias voltages showing the Coulomb staircase. 
Insets in (a) and (b) correspond to an array without charge 
or resistance disorder and respectively a = 0.5 and q = 1. 
The current shows a big jump at threshold and very weak 
small steps at higher voltages. The width of the steps is AE"' 
in (a) and 2EI"'' in (b). The staircase structure is completely 
washed out with increasing bias voltage. Main figures in (a) 
and (b) show I-V curves for arrays with charge disorder but 
homogeneous resistances, a = 0.5 in (a) and a = 1 in (b). 
The first step height is reduced and the staircase structure 
is more pronounced than in the clean case because tunneling 
processes in the bulk can have small energy gain and become 
the bottle-neck more easily. The width of the steps in (b) 
remains equal to 2EI"'' as in the clean case. In (a) there 
are two different step lengths which alternate and sum to 
iEc, corresponding to the entrance of charge through both 
electrodes. The I-V curves in (c) correspond to arrays with 
resistance disorder and a = 0.5 without charge disorder in the 
inset, and with charge disorder in the main figure. Disorder in 
resistance in general favors that bulk junctions act as bottle- 
neck for the current, but for a = 0.5 it affects differently to 
the steps associated to any of the contact junctions. 



in a bulk junction can sharpen the steps, as it creates a 
bottleneck for the current at a junction with an associ- 
ated energy for tunneling which does not directly depend 
on bias voltage, but the opposite behavior can also take 
place if the large resistance if found at any of the contact 
junctions. Note that the particular way in which the I-V 
curve is affected by disorder in resistance depends on the 
particular resistance and charge-disorder distributions. 



C. Linear Regime at High- Voltages 

At larger voltages linear dependence is recovered. The 
asymptotic linear I-V does not extrapolate to zero cur- 
rent at zero voltage, but it cuts the zero current axis at a 
finite offset voltage. This high-voltagc linear regime can 
be understood analytically. We obtain that the slope 
of the I-V curve depends only on the sum of the junc- 
tion resistances in series and the offset voltage is given 
by the sum of the excitonic energies of all the junctions. 
Contrary to what was found at low voltages, close to 
threshold, the offset voltage for purely onsite interac- 
tions recovers the zero inter-island capacitive coupling 
valueSS,!^ calculated starting from a finite value of the 
inter-island capacitance. To assist comparison with ex- 
periments we compute the voltage at which this linear 
regime is reached. 

At very high voltages, the charge gradient ensures 
that all the tunneling processes to the right decrease 
the energy. The corresponding tunneling rates are Ti = 
Ri^i<^t - Ef^''') and the total tunnehng rate r*°* for 
no resistance disorder is X^ta^Ti. Having in mind that 

ES'^f^ = r*°* = i?-' (^-ES'^r'')- This 

rate is independent on the selected tunneling process. To 
transfer a charge from the source to the drain requires in 
average {N + 1) tunneling events. The average current is 
thus 

4« (]vTT)i?^ " S ^^ J ^^^^ 

This prediction is compared in Fig. 4 (a) and (c) with nu- 
merical results. The slope of the current does not depend 
on a or the existence of charge disorder, but only on the 
number of junctions {N + 1). The slope is the same as 
obtained adding all the resistances in series. Asymptoti- 
cally, at low voltage, this curve cuts the 1 = axis at the 
offset voltage Vof/set = SS^-^i^"''- This value is, in 
general, different to the threshold voltage and indepen- 
dent of the resistance of the junction and the asymmetry 
of the bias voltage. 

Previous derivation is valid even if there is inhomo- 
geneity in the value of island capacitances but relies on 
the homogeneity of the junction resistances through the 
array. If this is not the case the total tunneling rate of 
each step in a sequence is r*°* = J2fJi^ ^»~^(*« - 
One could argue that on average, the charge gradient 
would be such that it ensures a uniform tunneling rate 
punj through all the junctions. The potential drop which 
gives such a tunnehng rate is = i?ir"''" + E^^'^ 

and F""* = E^''") /Rsurn with R,um = 

X^ila^ There are {N -f 1) possible tunneling events 
at each step in a sequence and (TV -I- 1) steps, thus both 
(A'^-l- 1) factors cancel our. The resulting average current 
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FIG. 4: (a) I-V curves calculated up to high-voltages for 
= 50, a = 0.5 for clean and charge-disordered arrays. Both 
curves approach the same asymptotic curve at high voltages, 
even if the threshold voltage is quite different. Theoreti- 
cal prediction is included for comparison. (b) Computed I-V 
curves for charge-disordered arrays with different junction re- 
sistances (solid lines) with their theoretical asymptotic pre- 
dictions (dashed lines). From top to bottom a curve corre- 
sponding to an array with all-equal junction resistances, an 
array with randomly assigned resistances and an array with 
all-equal random resistances, except the first one which is ten 
times Rt- The slope differs but all the curves have equal 
offset voltage, (c) Slope (main figure) and offset voltage (in- 
set) which give better fitting to the numerically computed 
current at high-voltages as a function of the number of is- 
lands in the array, for several array parameters, all of them 
with homogeneous junction resistances. For comparison the 
theoretical prediction (dashed-line) is included. The offset 
voltage is proportional to the number of islands and given by 
X^ilii^ while the slope goes like the inverse of the sum 

of the junction resistances added in series. It is inversely pro- 
portional to the number of junctions, when all resistances are 
equal. Independence on the value of a and the presence or 
absence of charge disorder is observed, (d) Voltage at which 
the high-voltage asymptotic behavior is reached, estimated as 
the value at which / — lasympt /I is smaller than a given value, 
1% and 5% in the figures.lt is slightly larger for disordered 
arrays, increases linearly with the number of islands and it is 
approximately three and 2.5 times the offset voltage. In large 
arrays a very large bias voltage can be required to reach this 
linear dependence. 



is 

lasympt ^ [ ^ ~ ^ 'j (1^) 



As in the uniform resistance case, the slope in the cur- 
rent corresponds to the addition in series of all the resis- 
tances. The predicted asymptotic high-voltage behavior 
is observed in Fig. 4 for arrays with different parame- 
ters. Note that as longer is the array larger voltages have 
to be applied to reach this voltage. The voltage Vimear 
at which the linear behavior is reached is estimated in 
Fig. 4(d). It is approximately three times the offset volt- 
age and slightly larger in the presence of charge disorder. 
In long arrays Vunear can become very large and the lin- 
ear behavior will not be easily reached experimentally. 



IV. POTENTIAL DROP THROUGH THE 
ARRAY 

The potential through the array can nowadays be 
measured^, but to our knowledge it has not been studied 
theoretically. In conventional ohmic systems with a lin- 
ear current-voltage relation V = IR the potential drops 
homogeneously through the array if the resistivity of the 
system is homogeneous. When the proportionality con- 
stant between voltage and current is given by the sum 
of the resistances in series, but these resistances are not 
all equal, the voltage drop at each point is proportional 
to the local resistance. The nanoparticle array I-V char- 
acteristics are highly non-linear and in general it is not 
obvious how the potential drops through it. At the is- 
lands the potential is the sum of the disorder and charge 
terms, while at the electrodes the potential is controlled 
by the applied bias. As seen above, there are two linear 
regimes which can be identified. At high voltages the 
differential conductance equals the inverse of the sum of 
the resistances in series. Naively, a potential drop pro- 
portional to the resistance at each junction could be ex- 
pected at these voltages, but we will see that this is not 
exactly the case. In the low-voltage regime the slope of 
the current does not correspond to the addition of the 
resistances in series, but it is proportional to a single one 
(to the sum of two of them in the symmetric case) and 
it is not clear that the potential drop should be propor- 
tional to the junction resistance. In this section we study 
the potential drop through the array at low, intermediate 
and high voltages and show that in none of these regimes 
the potential drop at a given junction is strictly propor- 
tional to its resistance. In the linear regimes we show 
that deviations from this proportionality are related to 
the offset voltages and filling of the array. In the inter- 
mediate regime the voltage drop oscillates with position 
and reflects correlations between the charges. 
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A. High- Voltage Regime 

We start with the high-voltage Unear regime, as it is 
the easiest to understand. Fig. 5 shows the average po- 
tential drop in a clean and a disordered array for a given 
bias voltage in the high-voltage linear regime. All the 
junction resistances are equal in the top figures. The av- 
erage voltage drop is equal in both the clean and disor- 
dered case, and at first sight it seems linear. A linear po- 
tential drop through the array implies a homogeneous av- 
erage junction potential drop$i. However, at the contact 
junctions $i is approximately E^^^ times smaller than at 
the bulk junctions. The voltage drop at each junction 
is not equal to the current divided by the junction resis- 
tance either, as could be naively expected. The reasons 
for these deviations can be found in the offset voltage in 
the TV curve and in Eq. ^ which gives the change in 
energy for tunneling. 

At high voltages the I-V curve is linear but the total 
voltage drop through the array does not equal Rsuml- 
As, seen in equation there is an offset voltage 

Voffset = T.t=i E^^^. This offset voltage reflects the 
excitonic energy cost for tunneling. The excitonic energy 
is not equal at each junction, ft is 2£'^*' at the bulk junc- 
tions and approximately E"^ at the contact ones. Only 
the extra potential drop $i — El~'^ at each junction gives 
a finite contribution to current through it. On average 



/ = 



1 



Rs 



- E. 



e—h\ 



(13) 



From current conservation at high voltages the average 
potential drop through the array 



= E' 



Ri 



R. 



{V -Voffset) 



(14) 



It is not affected by the presence of charge disorder in 
the array (but it would change if capacitances are not 
homogeneous, via E'r^'''). As observed in Fig. 5, Eq ^ 
gives a good estimate of the potential drop. The validity 
of Eq (fT4|) is better seen when $i — E'^~'^ is plotted. It is 
proportional to the resistance of the junction and equal 
in every junction if all resistances are the same. This 
statement is valid independently on the position of the 
resistance, as shown in the bottom figure of Fig. 5 and its 
inset and on the asymmetry a of the applied voltage (not 
shown). The dependence of $i on the junction resistance 
is easily understood. The tunneling probability through 
a junction depends on its resistance. It is inversely pro- 
portional to it. When the resistance is very large, the 
charge has a lesser tendency to jump from an island to 
its neighbor and it will spend more time in the island 
producing a dependence of the time-averaged potential 
drop on the junction resistance distribution. 

As seen above, in this high-voltage linear regime the 
current can be obtained from the average tunneling rate 
and correspondingly from the average potential drop. 
Deviations of the average value (5$i, i.e. the root mean 
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FIG. 5: Top: mainfigure shows the average potential at the 
islands as a function of position for a disordered array 
with N — 50, a = 0.5 and all junction resistances equal. 
The average potential drop through the array is close to be 
linear. As shown in the inset the average potential drop at 
the junctions $i is homogeneous only once the excitonic en- 
ergy is substracted. The value subtracted is smaller at the 
contact junctions where "1> is smaller. Error bars give an es- 
timation of the fluctuations of the potential drop. Bottom: 
main flgure (inset) show the average potential drop, with er- 
ror bars giving its root mean square, at the junctions with 
the excitonic energy substracted corresponding to a clean ar- 
ray with the first (middle) junction ten times larger than the 
rest. $i — -Ef"'' is proportional to the junction resistance. 
Note that this proportionality holds only once the excitonic 
energy is substracted. 



square (r.m.s.), are shown in the inset in Fig. 5, in the 
form of error bars. They are slightly smaller at the con- 
tact junctions as the potential at the electrodes is re- 
stored to its nominal value via a battery prior to any 
tunneling event and larger at junctions with a larger re- 
sistance. Fluctuations in the local voltage drop $i in- 
crease with applied bias voltage as the number of pos- 
sible charge states and the width of the distribution of 
hopping energies do. Fluctuations are larger at those 
junctions with a larger resistance, but (5$^/$^ is smaller. 



B. Low- voltage regime 

Close to threshold the current depends linearly on 
V — Vt- Here we show that in this linear regime the 
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average potential drop mainly reflects the charge state 
of the array at threshold. This charge state depends on 
the asymmetry of the voltage drop a and disorder and 
in a symmetrically biased clean array on the even or odd 
number of islands. For a = 1, charges enter from the 
source and (N — 1) bulk junctions prevent charge mo- 
tion. If an electron reaches the last nanoparticle, it can 
freely jump onto the drain at zero potential. There is 
no charge gradient at the drain junction. Consequently, 
the potential drop at this junction vanishes at threshold. 
On the contrary, at the (TV — 1) bulk junctions there is 
a charge gradient equal to unity, with the corresponding 
potential drop 2El''K To allow current Vb — </>! equals 
the excitonic energy of the first junction, approximately 
equal to El'^K Close to threshold, as the bottle-neck for 
the current is the entrance of electrons from the source 
the charge state of the array is most of the time equal 
to the one at threshold, and only perturbed by the fast 
passage of charges. The average potential drop, plotted 
in Fig. 6(a) for a clean array with all junction resistances 
equal is almost the same as the static potential drop at 
threshold. 

For a — 1/2 and odd number of particles, at threshold 
the charge gradient and charge potential drop in a clean 
array are respectively one and 2i?*'*' at the bulk junctions. 
At the contact junctions the potential drop is E^'^K As 
in the a = 1 case discussed above, the average potential 
drop in the linear regime close to threshold, is very close 
to the one found at threshold, which equals the excitonic 
energy at each junction, — E^^^ ~ 0, what can be 
seen in the main figure in Fig. 6(b). When the number 
of particles is even, the charge gradient at one of the 
junctions vanishes. As shown in the inset of Fig. 6(b), 

— E^^^ in the bottleneck regime is positive and equal 
for all the junctions. This reflects that every junction is 
uncharged with equal probability. 

In the disordered case, only those junctions with up- 
ward steps in the disorder potential are charged, and this 
is reflected in the average potential drop, in Fig. 6 which 
adds disorder, charge and bias potential. 

The threshold voltage does not depend on the resis- 
tance of the junctions, but the flow of charge does. This 
is reflected on the average potential drop on such very 
small scale that even if the threshold voltage potential 
drop is substracted at each junction and for reasonably 
large changes in resistance it is not visible (see main fig- 
ure in Fig. 6(d). This is different to the dependence 
observed in the high-voltage regime. For extremely large 
values of the resistance disorder a weak effect on the av- 
erage voltage at bias close to threshold can be seen (not 
shown). In this case, the potential drop at a junction with 
a larger resistance is slighltly larger than at the rest. At 
the adjacent junctions it is slightly smaller, what reflects 
the average charge state of the nanoparticles joined by 
the large resistance. 

The r.m.s of the junction potential drop at low bias 
voltages are very small in some of the cases analyzed, of 
order 10~^i?*^' in main figures in Fig. 6(a) and (b). Error 
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FIG. 6: Average potential drop tlirougli the array at bias volt- 
ages very close to threshold. Fluctuations are smaller than the 
symbols, (a) Clean N=50 array with a = 1. Notice that the 
potential drop vanishes at the last junction, as the source is at 
zero potential and there is no charge gradient at this junction. 
At the first junction is equal to the excitonic energy, which at 
a contact junction equals iS^^'. (b) Average potential drop for 
a symmetrically biased N=49 array in main figure (N=50 in 
the inset) with the excitonic energy substracted. Once the ex- 
citonic energy, associated to the charge gradient at threshold 
has been substracted the average potential drop almost van- 
ishes for N — 49. The homogeneous and positive value for 

= 50 reflects that every junction is uncharged with equal 
probability. (c)Average potential drop at the junctions $i 
corresponding to a disordered array with A = 50 islands and 
a — 1/2. (d) Main figure (inset) $ — E"^"'' corresponding to a 

= 50 clean array with the middle (first) junction resistance 
10.8 larger than the other ones. When the larger resistance is 
in the middle its efi'ect on the average potential drop at these 
small voltages is barely visible and the potential drop almost 
equals the one found in the inset in (b). This is the generic 
behavior found with resistance disorder very close to thresh- 
old. A special case is a N-even clean a = 0.5 array with the 
first resistance larger than the other ones, shown in the inset. 
The average potential drop differs considerably with respect 
to the one found in the inset in (b). The presence of the larger 
resistance at a contact junction modifies the average charging 
of the array and selects the opposite contact junction as the 
one which lacks charge gradient. 



bars are small because most of the time the array is at the 
threshold charge state. This is not the case in the inset in 
6(b) where fluctuations are or the order of O-SE^^**', or in 
the inset in 6(d) where they are small everywhere except 
at the last two junctions where they are of order E'^*'. 
Disorder also increases the fluctuations in the average 
voltages to values of the order of the charging energy. 
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FIG. 7: Average potential drop, with the excitonic energy 
substracted, #i — E^~'^ at the junctions as a function of po- 
sition at several values of the bias voltage, for which the cur- 
rent is in the Coulomb staircase regime, corresponding to a 
clean A'^ = 50 array and a = 1/2. From top to bottom 
V = 102, 104, 106, 108, 110, 112£;r' Curves have been verti- 
cally displaced to avoid overlap. The potential drop show 
almost regular oscillations which reflect a stationary (but 
not static) charge ordering, resembling the one expected in 
a Wigner crystal (see discussion in ref.-SS,). The number of 
maxima and minima does not change within a step and in- 
creases in two (for a = 1/2) from one steps to the next one. 
The position of maxima is slightly adjusted within a step. 
They start appearing close to the electrodes and reflect the 
entrance of charge from them. 



C. Intermediate Voltage Regime 

The most interesting regime to analyze the voltage 
drop is at intermediate voltages where the I-V curve show 
the Coulomb staircase. For the case of a clean array with 
capacitively coupled nanoparticles Stopa^ showed that 
the steps in the I-V characteristic correspond to alterna- 
tion of the charge density between distinct Wigner crys- 
talline phases. In our case, the interaction is short-range, 
but the possibility of such a Wigner state with charges pe- 
riodically ordered to minimize their repulsion, if present, 
could be observed via the potential drop. Periodic charge 
ordering should lead to oscillations in the potential drop 
along the array. Such an observation would be a clear 
evidence of correlated motion. 

The average potential drop (with the excitonic energy 
subtracted) through the array for several voltages corre- 
sponding to clean N=50 nanoparticle arrays is shown in 
Fig. 7. Clear oscillations are seen. Comparing the values 
of the bias voltage chosen with the position of the steps 
in the corresponding I-V curve in Fig. 7, the number of 
maxima/minima in the potential drop do not change in 
a given step in the Coulomb staircase. For symmetri- 
cally biased arrays, they increase in pairs from a step to 
the next one. For odd/even number of particles, there 
is always a minimum/maximum at the center of the ar- 
ray. The other maxima and minima tend to be as equally 
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FIG. 8: (a) and (b) show the average potential drop at the 
first and 25th junctions as a function of bias voltage for a 
clean array with A = 50 islands, a = 1/2 and no resistance 
disorder at intermediate voltages. The average potential drop 
at the contact junction shows clear oscillations as a function 
of the bias, while the potential drop at the center of the array 
depends monotonously on V. (c) and (d) show $i — E^~'^ 
for charge and resistance disordered 49-islands arrays respec- 
tively and Q = 1/2. The Coulomb staircase is much more 
pronounced in the presence of disorder. However, as shown 
here, the oscillatory potential drop structure found in Fig. 7 
and characteristic of Wigner-crystal like physics is destroyed. 



spaced as possible, but this is not exact. Inconmensura- 
bility beween the period of the oscillations and the lattice 
can distort equal spacement. Also, when new maxima or 
minima appear they are closer to the source and drain 
electrodes and move inwards, producing a movement of 
the other maxima and minima, with increasing voltage. 
This can be taken as a finite size effect of the Wigner crys- 
tal state. As the number of charges in the array increase 
with increasing bias voltage the amplitude and period of 
the oscillations decreases, approaching the high-voltage 
regime for which l>i — E'^^^ is homogeneous. 

The anomalous potential drop can be also seen in the 
potential drop at a given junction as a function of the bias 
voltage, shown in Fig 8(a) and 8(b) for junctions 1 and 25 
for a clean symmetrically biased 50-islands array At the 
first junction the potential drop show clear oscillations 
as a function of the bias voltage, which reflect that new 
charges state at the first island are allowed. The potential 
drop increases until an extra charge can be accumulated 
at the first nanoparticle, for larger voltages the average 
occupation of first island increases and the voltage de- 
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creases smoothly until a new value at which it increases 
again as the increase in occupation of first island cannot 
compensate the increase in the electrode potential. Os- 
cillations, but less regular and less pronounced due to the 
movement of maxima and minima discussed above, are 
also observed at intermediate junctions . Potential drop 
is much more homogeneous at the middle of the array, 
where there is always a minimum (or maximum) in the 
potential drop. 

Charge or resistance disorder alters the charge motion 
and frustrates the formation of this Wigner crystal like 
state, as seen in Fig. 8(c) and 8(d). This is the opposite 
behavior that would be naively expected if one just as- 
sociates the appearance of plateaux with the oscillations 
in voltage drop and emphasizes that the step profile is 
just a consequence of the dependence on the bias voltage 
of the tunneling rate of the processes which control the 
current. 

The r.m.s. of the junction potential drop is larger than 
at low voltages and smaller than at high voltages. It 
is of the order of the excitonic energy and reflects the 
variation in occupation of the island. It slowly increases 
with voltage. 



V. SUMMARY 

In this paper we have studied, analytical and numeri- 
cally, the current through an array of N metallic islands 
placed in between two large electrodes, the source and 
the drain at voltages Vb and Vat+i respectively. The ap- 
plied bias voltage V satisfies V = Vq — V^+i- Interac- 
tions are restricted to those charges in the same conduc- 
tor. The capacitive coupling between different conduc- 
tors vanishes. The nanoparticle level spacing is assumed 
negligible and transport is treated at the sequential tun- 
neling level. In the model used we take into account that 
the electrodes are not ideal voltage sources, but have a 
finite self-capacitance. This means that the voltage at 
the electrodes fluctuates in response to tunneling pro- 
cesses, but we assume that prior to the next tunneling 
event the nominal voltage is restored. Due to the large 
value of the electrode capacitance numerical results are 
barely modified and its effect is neglected in the discus- 
sion. The probability of each process depends on the 
change in energy involved, which can be written as the 
energy to create an electron-hole pair plus the difference 
in potential between the sites involved in the tunneling, 
see Eq.([3]). We have analyzed and clarified how the trans- 
port properties depend on the number of particles, the 
presence of charge or resistance disorder and the bias 
voltage, including how symmetrically it is applied. To 
quantify this symmetry we have introduced a parameter 
a as Vb = aV. 

We have shown that in the purely onsite interaction 
limit the dependence of the threshold voltage Vt of clean 
arrays on the number of nanoparticles differs qualita- 
tively of the dependence predicted for weakly coupled 



islands, as in the onsite case studied here a charge can- 
not flow freely through an empty array. For symmet- 
rically biased arrays Vr equals 2NE^^^ for odd N and 
2Ec{N —ly^ for even N . The even-odd effect disappears 
for forward biasing (a = 1,0) and Vt = Ei^'-{2N - 1), 
see Fig. 1(a). The threshold voltage is not affected by 
disorder in the junction resistances but it depends on 
the selected disorder configuration if charge disorder is 
present. With charge disorder, the average threshold 
voltage is independent on a and we recover previously 
predicted values Vt — Elf^N, see Fig. 1(b). 

At voltages marginally close to threshold, current is 
linear on (V — Vt) with a slope independent on the num- 
ber of particles N but which depends on the resistance 
of the contact junctions and on the bias asymmetry a, 
see Eqs. to ([9|) and Fig. 2. This dependence has been 
obtained both numerical and analytically and resolves 
previous controversy on the power-law close to thresh- 
old. It reflects that the junction through which charges 
enter into the array acts as a bottle neck. The range of 
voltages at which this linear dependence holds is proba- 
bly too small to be observed experimentally. Linearity is 
lost when the scattering rate for tunneling through the 
contact junctions is comparable to other tunneling pro- 
cesses scattering rates. 

The linear regime is followed by a Coulomb staircase 
at intermediate voltages. The width of the steps depends 
on a and on the presence of charge disorder. For clean ar- 
rays the bias voltage step width is 2i?J*' for forward bias 
and 4i?*^' for symmetric bias. The stepwidth changes if 
disorder is present but still depends on the value of a. 
The staircase profile depends on the junction resistances 
values. See Fig. 3. 

At high voltages current depends linearly on bias volt- 
age. The asymptotic I-V characteristic is given by 
and cuts the zero current axis at a finite offset voltage, 
see Fig. 4. The slope of the asymptotic linear I-V is given 
by the inverse of the sum of the junction resistances in 
series and the offset voltage takes into account charg- 
ing effects. The high-voltage linear behavior is reached 
for bias voltages Vimear approximately three times larger 
than the offset. Vimear can be very large long arrays. 

We have studied the potential drop through the array 
and showed that in none of the transport regimes stud- 
ied it drops completely linearly through the array. In 
the low-voltage linear regime the average potential drop 
mainly reflects the charge gradient created at thresh- 
old, see Fig. 6. The effect of disorder in resistances is 
extremely weak, except for symmetrically biased arrays 
with even number of particles. In the Coulomb stair- 
case regime the potential drop along the array show al- 
most periodic oscillations in disorder-free arrays which 
reflect the correlated motion and the tendency to form 
a Wigner-crystal like state, see Fig. 7. Such periodicity 
is destroyed by charge or resistance disorder in Fig. 8. 
In the high-voltage regime the ohmic dependence, and 
the associated proportionality between the potential drop 
and the junction resistance , is only recovered once the 
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excitonic energy is substracted, as seen in (|T4)) and Fig. 5. 
The mean value of the potential drop serves to compute 
the I-V characteristic in this regime. 

VI. APPENDIX: SIMULATION 

We numerically determine the time evolution of the 
state of an array of nanoparticles sandwiched between 
two large metallic leads symmetrically biased at poten- 
tials Vo and Vat+i, with Vq — Vat+i. The state of the 
array consists of the set of charges {Q/^} that occupy the 
array islands and the leads. The island charges take on 
integer values. The charge of the nanoparticles is modi- 
fied when an electron tunnels between two adjacent sites. 
The charges of the source and drain can take on any real 
value because they can be modified discretely via the tun- 
neling of charges or continuously via the charging of the 
leads by a battery. 

The evolution of the system is computed by means of 
a kinetic Monte Carlo simulation"*^. At each iteration 
a single tunneling event takes place. The time involved 
in this event t depends on the tunneling rates of all the 
possible tunneling processes. Each iteration starts from 
an initial charge configuration. First, it is computed the 
change in energy and the tunneling rate of the 2(A^ -t- 1) 
possible hopping events, corresponding to the tunneling 
of a single electron, to the left or to the right, through 
any of the (A^-|- 1) junctions. The probability of changing 
the initial configuration varies with time like 

pchange(^^^ ^ ^ _ pstay ^^^^ ^ ^ _ g-r*°*(t-to) (^^g) 

with to the time at which the preceding tunneling pro- 
cess took place and r*°* = Y^fJ^^^t and 
r~ are the tunneling rates through the i junction to the 
left or to the right, respectively, and are calculated from 
([1]). To sample the time interval between two hopping 
events we generate random numbers between [0, 1] to 
mimic P'^'^'^^se a,nd obtain t — t — to from ([TC)) . As the 
average of —lnP^^°'y is the unity, if one is interested only 
in the average values of the charge or the current, and 
not on its fluctuations, the time step could be fixed^^ to 
l/r*°*. This option is numerically faster. 

The relative probability of each tunneling event is 
r^/r*°'. To determine the hopping process which 



changes the charge state, the relative probabilities are 
consecutively arranged in the interval [0,1]. A second 
random number in this interval is generated to select the 
tunneling process. 

Then, the charge configuration is updated. After we 
modify the state of the system, we allow the external cir- 
cuit to return the leads to their applied bias values prior 
to the selection of the next hop. This effect is simulated 
by resetting the charges on the source and the drain to 
the values that restore the nominal applied bias. 

In order to remove all sensitivity to initial conditions, 
before we track the evolution of {Qi} as a function of time 
at any voltage, we perform N^q > 10* iterations to equi- 
librate the system. Following these iterations, we track 
the evolution of the charge state until the net number of 
electrons that arrive at the drain, Q drain, equals a very 
large number (> 10^). The average calculated current is 
given by 



ttot 

where ttot is the sum of all time intervals between hopping 
processes in the evolution runs. If a tunneling event in- 
volving the drain is selected, an amount = ±1 is added 
to Qdrain depending on whether an electron hopped to 
or from the drain. Current conservation ensures that the 
average current is the same through any junction. The 
minimum numbers of equilibration cycles, N^q, and evo- 
lution cycles (set by Qdrain) depend on the voltage. 

To calculate the average voltage drop, we assume that 
the system is in a given state a time equal to the interval 
T until the next tunneling event takes place. 
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